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We investigate the finite-temperature properties of an Anderson lattice with regularly depleted 
impurities. The physics of this model is ruled by two different magnetic exchange mechanisms: con¬ 
ventional Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction at weak hybridization strength V 
and a novel inverse indirect magnetic exchange (IIME) at strong V, both favoring a ferromagnetic 
ground state. The stability of ferromagnetic order against thermal fluctuations is systematically 
studied by static mean-field theory for an effective low-energy spin-only model emerging perturba- 
tively in the strong-coupling limit as well as by dynamical mean-field theory for the full model. The 
Curie temperature is found at a maximum for a half-filled conduction band and at intermediate 
hybridization strengths in the crossover regime between RKKY and IIME. 
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I. INTRODUCTION 


As has been pointed out by Nozieresf^^the presence of 
a correlated impurity in an a priori uncorrelated metal 
introduces effective interactions among the conduction 
electrons. The range of these interactions decisively de¬ 
pends on the strength of the impurity-host coupling. 
Consider the case of an Anderson impurity® 

H (4oa/cr + h.C.) 

{hi) <r=td cr=l,l 

+ c/(/;a-i/2)(/|4-i/2), (I) 

with annihilators Ci„,f„ referring to local conduction- 
electron and impurity orbitals, respectively. For the case 
of a Hubbard interaction U and a local hybridization 
V much stronger than the nearest-neighbor conduction- 
electron hopping t, an effective Hamiltonian with an 
almost local interaction characterizing the low-energy 
physics of the conduction-electron system can be de¬ 
rived explicitly.® This is achieved by means of degenerate 
fourth-order perturbation theory in the hopping terms 
which connect the neighbouring conduction-electron sites 
to the site io where the impurity is coupled to. To leading 
order, the effective model is given by 

ijAio 2 

HeS = -^ E E : (2) 
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where z is the coordination number of the lattice, where 


a = t 
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is the effective interaction strength, and where S'bond is 
the spin-operator referring to the “bonding” symmetric 
linear superposition of the z orbitals neighbouring zq (see 
Ref.for details). 


There are three different energy scales to be consid¬ 
ered: (i) Local singlet formation at io takes place on 
the high-energy scale ~ U,V. While this singlet may 
be called a local Kondo singlet, its binding energy scales 
linearly with V for strong V. This is opposed to the 
weak-coupling limit 1/ —>■ 0 (with U t fixed) where it 
is exponentially small and where t he l ow-energy physics 
is dictated by a single Kondo scale.®^ (ii) On an energy 
scale ~ t, conduction electrons scatter at the local Kondo 
singlet. This scattering effect is already included at ze¬ 
roth order in the perturbative expansion and is formally 
described by excluding the site Zq from the summation 
in the first term of the effective Hamiltonian in Eq. Q. 
(iii) The first non-trivial effect takes place at fourth order. 
An effective interaction among the conduction electrons 
in the immediate vicinity of the impurity emerges which 
is mediated by virtual excitations of the local Kondo sin¬ 
glet. This happens on the lowest energy scale given by 
the effective coupling constant a in the second term of 
Eq. ([2|. 

A fundamentally interesting question is whether the 
emergent effective interaction among the a priori uncor¬ 
related conduction electrons can give rise to collective 
phenomena. This may be expected for a lattice variant 
of the model, i.e., for a system with a thermodynamically 
relevant concentration of impurities. The extreme case 
is a periodic Anderson model with a depleted system of 
“impurities” placed at every second site, i.e., on the B 
sites of a bipartite lattice consisting of sublattices A and 
B. Fig. [l] displays an example for the D = 3 dimensional 
simple-cubic lattice. We consider a model with L sites 
(L —> oo in the thermodynamical limit) and R = L/2 
impurities. The total number of electrons N satisfies 
2R < N < AR such that there are well-formed local 
Kondo singlets in the low-energy sector. 

At fourth order, perturbation theory is essentially un¬ 
changed as compared to the impurity model Eq. Q, since 
any local Kondo singlet, consisting of the correlated im¬ 
purity coupled to an B-sublattice site, is surrounded by 
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uncorrelated A-sublattice sites, and thus the same virtual 
processes lead to the same effective interaction. There¬ 
fore, the resulting effective Hamiltonian only involves A- 
sublattice sites and excess conduction electrons that are 
not absorbed in a local Kondo singlet. The hopping term 
becomes ineffective since the excess conduction electrons 
are confined between the local Kondo singlets surround¬ 
ing each A site. Hence, we are left with a lattice model 
of A sites, made up by non-local spins Si^hond referring 
to the bonding orbital around each B site: 

-^eff = ^ ^ '^i>ond J (4) 

ieB 

with 5i,bond = where cr is the 

vector of Pauli matrices and where the creation operator 
of the bonding orbital around i G B is given by = 
i.e., the bonding one-particle orbital is 
the symmetric superposition of neighboring A-sublattice 
orbitals: 

n.n.{i) 

Ibondi,cr)= ^ |j,cr)/v^ (iGB). (5) 

j&A 

The effective spin-only model, Eq. 0 , is non-trivial as 
different non-local spins refer to overlapping orbitals and 
therefore do not commute. 

There is not much known about this model: At half¬ 
filling, N = L + R = 3i?, one can rigorously show that 
a ferromagneticP state with fully polarized magnetic mo¬ 
ments of the conduction electrons on the A sites is among 
the ground states.!^ Exact diagonalization of small sys¬ 
tems suggest^ that the model has a ferromagnetically 
ordered ground state in the filling range 2R < N < AR 
(for lower or higher fillings, local Kondo singlets are bro¬ 
ken up). An inverse indirect magnetic exchange (HME) 
where the magnetic moments of A-site electrons are cou¬ 
pled ferromagnetically via virtual excitations of the local 
Kondo sin glets has been identified as the main physical 
mechanism.l21i21 For a one-dimensional depleted Ander¬ 
son lattice, density-matrix renormalization-group calcu¬ 
lations have showrP that the HME mechanism gradu¬ 
ally crosses over to a conventional (RKKY)PiE^ indirect 
magnetic exchange, also favoring ferromagnetism, when 
varying V from strong to weak hybridization at fixed 
U t. This crossover and the mutual interplay between 
RKKY and HME mechanisms for the magnetic ground- 
state properties has recently been discussed in Ref. |T3] in 
the context of SU(N) models of ultracold Fermi atoms 
trapped in optical lattices. 

The purpose of the present paper is to study the finite- 
temperature properties of the depleted Anderson lat¬ 
tice, particularly the stability of the ferromagnetic or¬ 
der against thermal fluctuations. From the RKKY the¬ 
ory, one can expect Tc (x J'^ oe for the Curie tem¬ 
perature at weak V and in a pa rame ter regime where 
the Schrieffer-Wolf transformatiorP^li^ applies such that 



FIG. 1: (Color online) Depleted periodic Anderson model 
with R = L/2 impurities on a D = 3-dimensional simple- 
cubic lattice with L sites {L -A oo). Correlated impurity sites 
(red) with on-site Hubbard interaction U are coupled via a 
hybridization of strength V to the B sites (blue) of the bi¬ 
partite lattice. For strong U,V 3> t, local Kondo singlets are 
formed on the half-filled “dimers” consisting of impurity and 
B sites (if the total electron number N satisfies 2R < N < AR) 
and strongly conhne the motion of the excess conduction elec¬ 
trons on the A-sublattice sites (green). Virtual excitations of 
the local Kondo singlets induce an effective interaction of the 
conduction electrons on A sites. 


J = jU. On the other hand, for strong V, the ef¬ 
fective model Eq. Q suggests that Tc oc a oc We 

therefore expect a pronounced maximum of Tc at an in¬ 
termediate V. This optimal V but also the absolute value 
of Tc are interesting from a fundamental theoretical per¬ 
spective. Not only the strong V dependencies but also 
the fact that the non-interacting {U = 0) depleted An¬ 
derson lattice exhibits a flat band at the Fermi energjffSl 
promise a comparatively high value for the critical tem¬ 
perature. Furthermore, the finite-temperature properties 
are important for the question whether magnetic corre¬ 
lations and magnetic long-range order induced by the 
HME can be verified experimentally. Candidate systems 
are magnetic nanostructures on non-magnetic surfaces as 
their geometrical and magnetic properties can be mea¬ 
sured, controlled and manipulated to a high degree on a 
atomic scale.^i^t^ Likewise, ultracold-atom systems come 
into question, due to the rapidly improving experimental 
techniques in this field and particularly due to the recent 
advances to employ fermionic alkaline-earth atoms to ef¬ 
ficiently simulate systems with spin and orbital degrees 
of freedom.l22H2ll 

Our study is based on two different types of mean-field 
methods: To address the strong-K limit, we apply static 
mean-field theory to the effective spin model Eq. 0. 
Since S'^^bond is not a rigid spin with S = 1/2, a fermion 
mean-field approach must be employed. Using this ap- 
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proximation, a rough estimate of the dependence of the 
Curie temperature on lattice dimension or coordination 
number and electron density is obtained. S econd ly, we 
apply dynamical mean-field theory (DMFTj^MSI to the 
depleted Anderson lattice. For a model with a de¬ 
pleted system of correlated sites, the DMFT can ex¬ 
pected to yield reliable results since the electron self¬ 
energy is much more local as compared to the dense 
model. This has been checked for the D = 1 dimen¬ 
sional model where essentially exact results are available 
via the density-matrix renormalization group technique.^ 
For ground-state properties of local observables as ob¬ 
tained by DMFT even quantitative agreement has been 
found. 

The paper is organized as follows: The static and dy¬ 
namical mean-field methods are introduced along with a 
discussion of the corresponding results in Secs.|n]and III 
respectively. The conclusions are summarized in Sec. IV 


II. STATIC MEAN-FIELD THEORY 
A. Depleted Anderson lattice 

The Hamiltonian of the depleted Anderson lattice is 
given by 

c\aCia + V -h h.C.) 

+ + (e - m) Y ■ 

i,cr i£B,(T 

( 6 ) 

It describes a system of electrons hopping over the sites 
of a bipartite D-dimensional lattice consisting of L sites 
with periodic boundary conditions. The two sublattices 
are denoted by A and B. We consider a D = 3 simple- 
cubic lattice (see Fig. but also the corresponding one- 
and two-dimensional cases (Figs. and [^. creates 
a conduction electron in a one-particle orbital with spin 
projection a at the site i = The nearest- 

neighbor hopping t = 1 sets the energy scale. 

One-particle orbitals at the B sites of the lattice hy¬ 
bridize with orbitals at R = Lj^ additional “impurity” 
sites with hybridization strength V . creates an elec¬ 
tron at the impurity site attached to site i G B oi the sub¬ 
lattice B. Furthermore, and n’'/J = flf^^ 

denote the occupation-number operators for A, B and 
for impurity sites, respectively. The / orbitals should be 
considered as magnetic orbitals: There is a finite repul¬ 
sive Hubbard interaction U and the one-particle energy 
is set to e = —U/2 such that, for strong U, the forma¬ 
tion of local magnetic moments at the impurity sites is 
favored. 

The Hamiltonian Eq. (§ contains an overall chemical 
potential /r, i.e., we work with the grand canonical ensem¬ 
ble where /r is used to fix the average number of particles 


I . I . I . 

A B 

FIG. 2: (Color online) Geometry of the one-dimensional 
dilnted Anderson lattice. Red: “Imparities” with finite Hub¬ 
bard interaction. Green and blue: sites of the A and of the B 
sublattice, respectively. 



FIG. 3: (Color online) The same as in Figs. [^and Fig. for 
the two-dimensional case. 

{N). We will consider the range 2R < (N) < AR for our 
calculations. 

Switching off the hopping, i.e. t = 0, defines an atomic 
limit of the model Eq. (§. The ground state in the 
atomic limit is highly degenerate. For the considered 
range of the total electron number, each ground state is 
characterized by completely local Kondo singlets formed 
on the B and the attached impurity sites binding two elec¬ 
trons per singlet. The ground state degeneracy is due to 
the various configurations of remaining electrons on the 
A sites. Their density ua = /La, where 

La = L/2 \s the number of A sites, can vary within the 
range 0 < riA < 2. 

The depleted Anderson lattice Eq. (§ exhibits the con¬ 
ventional U(l) and SU(2) symmetries corresponding to 
conservation of the total particle number and the total 
spin. For /r = 0 the system is half filled, i.e., {N) = 3i? 
or riA = I, and there is an additional SU(2) isospin 
symmetry.l^ Due to particle-hole symmetry, we can re¬ 
strict our considerations to the range at and below half¬ 
filling. 

B. Strong-coupling limit 

For strong V ^ t, an effective Hamiltonian ’Heg can 
be derived by means of fourth-order perturbation theory 
in t around the degenerate atomic limit .1^ In this limit 
the ground state is characterized by local Kondo singlets 
at the B sites and a residual low-energy dynamics of the 
A-site electrons which is mediated by virtual high-energy 
excitations of the local Kondo singlets. Hence, HeS con¬ 
tains A-site degrees of freedom only. There is a very com¬ 
pact and highly symmetric representation of HeS given 
by Eq. @ with the coupling constant a specified by Eq. 
([^. Details of the perturbation theory can be found in 
Ref. [SI 

Here, we rewrite the effective Hamiltonian such that a 
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static mean-field decoupling can be applied in a straight¬ 
forward way. To this end, we use the dehnitions given 
below Eq. 0 to express the non-local spin operators in 
terms of creators and annihilators for electrons on A sites. 
Furthermore, we switch to a representation in reciprocal 
space by means of Fourier transformation in the form 

= V (zeA) (7) 

^ kGBZA 

where La = L/2 and where fc is a wave vector in the 
Brillouin zone BZa of the reciprocal A sublattice. Note 
that the A sublattice is a square lattice for L) = 2 but a 
b.c.c. lattice for the D = 3 case with a unit cell spanned 
by the basis vectors 01 , 02,03 displayed in Fig. With 
this we get: 

'HeS = ^^(A'(k) — ^)cj^^Ck(T 

k 

+ ^ ^PqkCp.^Cp_ktCq^Cq+k;^ . (8) 

p,q,k 


mean-held Hamiltonian which allows for a straightfor¬ 
ward study of the temperature dependence of the A-site 
magnetic moment and therewith gives access to the criti¬ 
cal (Curie) temperature Tq. However, typical mean-held 
artifacts must be expected and tolerated. 

Using the decoupling 0 in Eq. ([^ , we obtain a mean- 
held Hamiltonian 

'HeS = ( 12 ) 

k,(T 


which is bilinear in , c. The mean-held dispersion 


ry.(k) = --(U-Q_.) 7 ^(k) 


(13) 


as well as the constant in Eq. (12) depend on the possibly 
spin-dependent mean held which must be determined 
self-consistently from the following mean-held equation: 




1 


g/3(r;„(k)-/i) _|_ I 


(14) 


The effective one-particle dispersion is given by: 

T)f\ 

£;(k) = -—7^(k) (9) 

where £o(fc) = —^{k)t is the tight-binding dispersion of 
the L>-dimensional lattice. This also determines the k 
dependence of the interaction parameters of the effective 
Hamiltonian via: 

Up,q,k = |7(p)7(q)7(P - k) 7 (k + q) . (10) 

Apparently, the effective Hamiltonian describes itinerant 
electrons on the A sublattice with an interaction, the k 
dependence of which corresponds to the non-locality of 
the quartic parts of the Hamiltonian in real-space repre¬ 
sentation Eq. Q. 

C. Mean-fleld approximation 


Here /3 = 1/T and we have chosen units such that fcs = 1. 
The spin-dependent average A-site occupation number, 
n-Acr = I i® obtained as 


_ J_ \ _ 

~ La ^ -I-1 

k 


(15) 


With this, the order parameter, i.e., the A-sublattice 
magnetization, is given by toa = "^At ~ 

Numerical calculations are performed by starting with 
a guess for the chemical potential and solving the coupled 
system of Eqs. (131 and (141 self-consistently for each 
spin projection. From the self-consistent mean held Qo-, 
we obtain riAa via Eq. (15). In an outer self-consistency 
loop we then adjust the chemical potential until the to¬ 
tal hlling nA\ + u-a^ equals the given hlling tla- In the 
case of half-hlling ua = 1 , calculations are facilitated by 
particle-hole symmetry which hxes the chemical potential 
to /r = 0 . 


Note that in the strong-coupling limit both, the one- 
particle part as well as the interaction, scale with a. 
Therefore, the standard mean-held decoupling of the in¬ 
teraction term, 

'^pt^P-kt^q^'^k+q^ —>■ (Cp^Cp_kt)Cq^Ck-(.q^ 

+ Cp.^Cp_kt(Cq^Ck-l-q^) 

— (Cp^Cp_kt)(Cq^Ck-|.q^) , (H) 

cannot be controlled by a small parameter but must 
rather be seen as a Hartree-Fock approach neglecting 
correlation effects in the low-energy sector and assuming 
a collinear and homogeneous structure of the magnetic 
moments. The formal advantage is that one obtains a 


D. Results 

Calculations have been performed for lattices with dif¬ 
ferent dimensions D — 1,2,3 (see Figs. mil and re¬ 
spectively) as well as for different hllings ua at and be¬ 
low half-hlling. Fig. shows the resulting self-consistent 
mean helds Qa as functions of the temperature. For any 
D and riA, there is a non-zero critical temperature Tq 
below which we hnd a spontaneous spin-splitting of the 
mean held. This supports the above-mentioned exact- 
diagonalization results of Ref. 0 where a fully polarized 
magnetic ground state has been found for small one¬ 
dimensional systems in the hlling range considered. 

As can be seen in Fig. there is only a weak de¬ 
pendence of the mean held on the dimension D - after 
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FIG. 4: (Color online) Spin-dependent mean field [see Eq. 
(14l] as a function of the reduced temperature T/Tc for dif¬ 
ferent dimensions D (top, middle and bottom panel) and dif¬ 
ferent fillings ua as indicated (see top panel) at and below 
half-filling (tia = 1). Solid lines: cr Dashed lines: a =4,. 


rescaling with D or with the coordination number 
z. For r = 0, we have Qi = 0 and thus the a 
mean-field dispersion simplifies to = —aD^‘^{k)/2 
resulting in = max. and, at half-filling, = z 
since particle-hole symmetry enforces — D = D — Q^. 
For T higher than the Curie temperature Tc, we have 
Qt = Qi- The spin-independent mean-field is slightly 
decreasing with increasing T, except for half-filling where 
= Qi = zl2 = const, above Tq- 

Fig. [5] shows the temperature-dependent magnetiza¬ 
tion for the different fillings and dimensions. At zero 
temperature, the system is always fully polarized, i.e., 
nAi = 0 and riAt = Similar to the discussion of 
the mean fields, after rescaling the temperature with the 
respective Curie temperature Tc, there is a weak depen¬ 
dence of TOA on the dimension D at finite T. The phase 
transition from the ferro- to the paramagnetic state at 
Tc is of second order for any ua- Close to the Curie 
point, we find a critical behavior of toa characterized by 
the (mean-field) critical exponent for the magnetization 
P = 0.5, as expected. 

The only unexpected result consists in the unconven¬ 
tional T-dependence of toa at half-filling. While at low 
temperatures the missing feedback of long-wavelength 
spin excitations explains the absence of a power-law 
T dependence, one would expect, as a typical mean- 
field behavior, an exponential convergence of toa(T) for 


FIG. 5: (Golor online) Order parameter toa = n-At — u-aj. [see 
Eq. ( |15| )] as a function of T/Tc for different dimensions D 
and fillings ua- 


T —^ 0 with a negative curvature and a vanishing slope 
limr_,,o dTOA/dT = 0. However, for ua = 1, Fig. [^shows 
an inflection point of toa(T) at a finite temperature, 
which is increasing with increasing D, and an unusual 
upturn of Toa for T —>■ 0. Closer inspection of the data 
shows that the slope is diverging: 


dniA 1 

dT ^ y/odlT 


and 


dm A ln(T/a) 
dT ^ y/PdjT 


(D = = 3) 


(T = 2). 


(16) 

(17) 


The reason of this behavior is a van Hove singularity of 
the spin-dependent mean-field local density of states at 
the Fermi edge and is discussed in Appendix fA} 

From the temperature-dependence of the order param¬ 
eter we can read off the Curie temperature. This is plot¬ 
ted in Fig. for different D as functions of the filling ua- 
Since Tc(nA) = Tc(2 — ua) due to particle-hole symme¬ 
try, we restrict ourselves to the range ua < 1. Clearly, 
the Curie temperature must be proportional to the cou¬ 
pling constant a as there is a single energy scale in the 
effective Hamiltonian Eq. Q. 

Its geometry dependence is more interesting: Namely, 
Tc is by no means proportional to the coordination num¬ 
ber as it typical for many mean-field approaches but is 
much more rapidly increasing with increasing z (note 
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FIG. 6: (Color online) Filling dependence of the Curie tem¬ 
perature for lattices with different dimensions as obtained 
from the static mean-field theory. Note that Tc is rescaled 
by and given in units of the coupling constant a. 


that the numerical results are scaled by a factor in 
Fig. |§. This finding is not related to singularities in 
the density of states as it holds for any filling. We at¬ 
tribute the unconventionally high Tq to the non-locality 
of the interaction in the effective Hamiltonian and to the 
resulting fe-dependent contribution of the mean field to 
the mean-field dispersion in Eq. (13). 


III. DYNAMICAL MEAN-FIELD THEORY 
A. General theory 

Dynamical mean-field theory (DMFT^PSEfil neglects 
the feedback of non-local, e.g. magnetic correlations, on 
the local self-energy and the local one-particle Green’s 
function but correctly accounts for all local correlations. 
This represents a decisive step beyond the static mean- 
field approach. Particularly, the DMFT is able to de¬ 
scribe the formation of local magnetic moments already 
in the paramagnetic phase of a lattice model of itiner¬ 
ant electrons, such that the phase transition between the 
paramagnetic and the ferromagnetic phase at Tc can be 
understood as a transition between well-formed but dis¬ 
ordered moments and long-range order. This is opposed 
to the static theory where the local moments essentially 
vanish above the Curie point. 

It is important to note that the feedback of non-local 
correlations neglected within single-site DMFT is much 
weaker for the depleted Anderson lattice considered here 
as compared to a lattice fermion model with a dense sys¬ 
tem of correlated sites. This can be understood in the 
following way: Formally, the only approximation to be 
tolerated within DMFT is the locality of the self-energy. 
For a dense lattice model, such as the Hubbard model, 
for example, the DMFT becomes exact in the limit of 
infinite spatial dimensions D ^ oo since the nearest- 
neighbor elements of the self-energy scale as as 


can be inferred from its diagram expansion.!^ This is re¬ 
lated to the scaling of the bare propagator, 

where d = ||f — j|| is the Manhattan distance between 
the orbitals at sites i and j of a hyper-cubic lattice of 
dimension D. For a depleted Anderson lattice on a high¬ 
dimensional bipartite lattice with a Manhattan distance 
d between the correlated sites, this also implies that the 
non-local elements of the self-energy exponentially dimin¬ 
ish with increasing d. A completely local self-energy is 
realized in the single-impurity limit d —)■ oo. For the Hub¬ 
bard model (d = 1) and the periodic Anderson model 
(d = 3) and for low dimensions, quantitative studies 
have been performed within second-order perturbation 
theory.l^aHaSl gg^gg studied here corresponds to d = 4 
but there are two, possibly largely different hopping pa¬ 
rameters, t and V. For the ground state of the depleted 
Anderson lattice (with d = 4) in D = 1 dimension, a 
direct comparison between DMFT and essentially exact 
results obtained by the density-matrix renormalization 
group (DMRG) method has been performed in Ref. [HI 
and excellent agreement has been found for static local 
observables in the entire V/t regime. Comparing with 
DMRG, a quantitative discussion of the artifacts of the 
DMFT has been given in Ref. for a D = 1 tight- 
binding model with two Anderson impurities. Conclud¬ 
ing, we therefore expect that the DMFT yields reliable 
results. 

DMFT is easily adapted to the model Eq. ^ : For any 
dimension D, there are three sites in a primitive unit cell 
of the lattice (see Figs. BE and . Hence, the single¬ 
particle Green’s function Gk{uj) is a 3 x 3-matrix for any 
wave vector k in the first Brillouin zone of the A sublat¬ 
tice and for any one-particle excitation energy w. Sum¬ 
mation over k provides us with the local Green’s func¬ 
tion with, say, the (3,3) element referring to the impurity 
Green’s function {{fia', fja))uj- Using Dyson’s equation, 
this can be obtained from the local self-energy E(w) as 


G 


loc 


feeBZ^ 


uj + fi — e{k) — S(a;) 


(18) 


oil3 


Here, a,/3 = 1,2,3 label the different sites in a unit 
cell. Furthermore, S(a;) is a 3 x 3 diagonal matrix with 
5133 ( 0 ;) = E(o;) and Sii(o;) = £ 22 ( 0 ;) = 0 , and 

/ 0 eo(fc) 0 \ 

e{k) = eo(fc) 0 U (19) 

\ 0 V £ ) 


is the lattice Fourier transform of the hopping parameters 
with eo(fc) = —'y{k)t. 

The DMFT self-energy E(o;) is obtained as the impu¬ 
rity self-energy of an effective Anderson impurity model 
specified by the Hubbard-f7 and a hybridization function 
that is fixed by the self-consistency equation of DMFT 
as 

A{uj) = uj + n - £ - T{uj) - I . (20) 

Gloc (w) 
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Here, the impurity one-particle energy is given by e, and 
S(a;) must be determined self-consistently with Eq. (181. 


To compute the self-energy of the effective impu¬ 
rity problem at finite temperature T, we e mploy the 
continuous-time quantum Monte-Carlo methocP^^^Sl 
the hybridization expansion of the action of the effective 
impurity modelP^ Configurations are samp led by means 
of the Metropolis-Hastings algorithm PSM] As the Hub¬ 
bard interaction is of density-density type, we can use the 
highly efficient segment-picture variant of the approach 
and, following Ref. HU directly measure the impurity self¬ 
energy Tia(iujn) on the fermionic Matsubara frequencies 


B. Results 
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From the results of static mean-field theory for the ef¬ 
fective low-energy model Eq. (|^ in the strong V limit 
(see Fig. we infer that the Curie temperature is at a 
maximum for half-filling. More generally, we expect that 
at half-filling the stability of a ferromagnetically ordered 
state against thermal fluctuations is the highest not only 
for strong V but also for weak V, where the period of 
the RKKY interaction is commensurate with the posi¬ 
tions of the correlated sites on the lattice. Furthermore, 
at half-filing and for the considered lattice geometries, 
the RKKY interaction is ferromagnetic. We will there¬ 
fore restrict ourselves to the particle-hole symmetric case 
with the chemical potential fixed at /x = 0 and with the 
one-particle energy of the impurities set to £ = —{7/2 [see 
Eq.|6])]. We also fix the Hubbard interaction at an inter¬ 
mediate value {7 = 8 for the rest of the paper. To discuss 
the crossover from the RKKY limit to the regime of the 
inverse indirect magnetic exchange, we consider different 
hybridization strengths V. 

By carrying out a sum over Matsubara frequencies, one 
may easily compute the average spin-dependent occupa¬ 
tion numbers on the A and B sites, UAa = (cl^^CAcr} and 
n-Bcr = (cg^CBcr), from the local Green’s function Eq. 
(d), once self-consistency has been achieved. The av¬ 
erage occupation numbers of the impurity site, riimp o- = 
iflU). can be obtained in the same way or, equivalently, 
can be measured within CT-QMC directly. 

We explicitly allow for symmetry-broken states with 
finite ordered magnetic moments toa, xub and Wimp, 
where we have defined uia = riA^ — nAi, and ttib, mimp 
analogously. It is found that magnetic solutions of the 
DMFT equations are easily stabilized in the entire range 
of hybridization strengths V considered (but for suffi¬ 
ciently low temperatures) by starting the self-consistency 
cycle with a slightly spin-asymmetric initial self-energy. 
In addition, we also compute the homogeneous static im¬ 
purity spin susceptibility of the paramagnetic state x = 
5”T-imp/9i3|B=o- Here, B is the strength of a homoge¬ 
neous magnetic field coupling to the z-component of the 
total impurity spin as H H — B ifltM - flfa) 
where H is given by Eq. ([^. 


FIG. 7: (Color online) Ordered magnetic moments jtia, me, 
uiinip (circles) on the A sites, the B sites and the impurity 
sites, respectively, and the inverse homogeneous static im¬ 
purity magnetic susceptibility x~^ (diamonds) as functions 
of temperature T as obtained by DMFT for the D — 3- 
dimensional depleted Anderson lattice (see Fig. [^. Hubbard 
interactiou: {7 = 8, hybridization strength: V = 2. The line 
indicates a linear fit to the trend of x~^{T). The temperature 
and energy scales are fixed by the nearest-neighbor hopping 
t = 1 [see Eq. ([^]. 


Fig. 0 shows the results of a DMFT calculation at 
17 = 2 for the D = 3-dimensional lattice with L = 52^ 
sites, with additional R — L/2 impurites and periodic 
boundary conditions (see Fig. [^. This is fully sufficient 
to ensure that the results do not significantly depend on 
L. Statistical errors of the quantities shown in this and in 
the following figures are smaller than the size of the sym¬ 
bols. A typical Monte-Carlo run consists of more than 
10^ sweeps, and each sweep of more than k Monte-Carlo 
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FIG. 8: (Color online) The same as in Fig. but for V = 3. 












FIG. 9: (Color online) rriA, ms, mimp (circles) and x~^ (di¬ 
amonds) as functions of T, as in Fig. but for the D — 1- 
dimensional depleted Anderson lattice ^ee Fig. |^. Hubbard 
interaction: U = 8, hybridization strength: V = 2. DMRG 
data for T — 0 (squares) from Ref. [9] are included for compar¬ 
ison. 


steps with k being the average expansion order. Less 
than 50 DMFT iterations are sufficient for convergence 
of the results within the statistical error. 

For high temperatures the system is in a paramag¬ 
netic state. The inverse susceptibility x~^ shows a lin¬ 
ear Curie-Weiss trend from which one can safely esti¬ 
mate the value for the Curie temperature Tc « 0.064. 
X is calculated from the magnetic moments induced by 
an explicitly applied homogeneous field for sufficiently 
weak field strengths in the linear-response regime (typi¬ 
cally B < 0.01). 

The transition to the ferromagnetic state at low tem¬ 
peratures appears to be of second order, and the data for 
the ordered magnetic moments are consistent with a lin¬ 
ear temperature trend of m? close to Tc, i.e. rri^ cx (Tc — 
T). This implies a critical exponent /3 = 0.5 as it must 
be expected for a DMFT calculation. Note, however, 
that due to critical slowing down, it becomes progres¬ 
sively more difficult to stabilize symmetry-broken DMFT 
solutions for temperatures close to Tc- The double oc¬ 
cupancy at the impurity site dimp = ('raimp,t^imp,^), and 
thus the local magnetic moment = 3(1 — 2di„ip)/4 
turns out to be almost constant in the entire tempera¬ 
ture range considered: dimp ~ 0.077. In particular, the 
moment does not change signihcantly across the phase 
transition. The transition to the symmetry-broken state 
must therefore be characterized as long-range ordering 
of local magnetic moments that are preformed at higher 
temperatures. This is a typical effect of strong correla¬ 
tions and opposed to simple Hartree-Fock-like (or Stoner- 
like) phase transitions where the local magnetic moment 
forms right at Tq. 

The low-temperature state of the system actually dis- 


FIG. 10: (Golor online) The same as in Fig. but for V = 3. 

plays ferri-magnetic order since the magnetic moment 
at the B sites is antiferromagnetically aligned (tob < 
0) to the moments at the impurities and the A sites 
(tninip,mA > 0). This alignment is reminiscent of the 
antiferromagnetic coupling in the Kondo limit of the 
model, i.e. for R —>■ 0, where an antiferromagnetic effec¬ 
tive exchange interaction (Kondo coupling) of strength 
J = SV^jU emerges betwe en B sites and impurities 
in the low-energy sector.l^^li^ In the weak-coupling limit 
K —>■ 0, one furthermore expects that well-formed lo¬ 
cal magnetic moments appear at the impurity sites since 
charge fluctuations are strongly suppressed. Ferromag¬ 
netic coupling of these moments via the RKKY exchange 
then implies Imimpl —>■ I, while mA,mB 0. For 1/ = 2, 
we are still in the RKKY regime since the A-site moment 
is clearly smaller than the moment on an impurity site. 

As Fig. demonstrates, however, this changes with 
increasing V. For V = 3, we find uia > niimp at low 
temperatures indicating the crossover from the RKKY 
regime to the strong-K limit. In the strong-coupling limit 
V ^ t, almost localized “Anderson singlets” are formed 
by the magnetic moments at B and impurity sites, and 
thus wb, minip —5^ 0. The presence of local singlets at the 
B sites implies that electrons on the remaining A sites 
are very efficiently localized such that well-formed local 
moments emerge. Those moments couple ferro magn eti- 
cally via the inverse indirect magnetic exchange , “21121 j 
by virtual excitations of the Anderson singlets, and thus 
niA —>■ 1. This picture well explains that iriA > Wimp in 
Fig. I 

It is instructive to compare the results for the D = 3 
lattice with those obtained for D = 1 (see Fig. [^. Figs. 

and [T0| show results for the ordered magnetic moments 
and the impurity magnetic susceptibility for a chain ge¬ 
ometry with L = 50 sites with periodic boundary condi¬ 
tions. The overall trends seen in the figures are similar 
to those found for D = 3 but the crossover from the 
RKKY to the IIME regime appears at lower hybridiza- 
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FIG. 11: (Color online) Cnrie temperature Tc for the D = 3- 
dimensional depleted Anderson lattice at 1/ = 8 and half¬ 
filling as a function of the hybridization strength V. Points are 
obtained via X~^{Tc) = 0 by extrapolating the linear temper¬ 
ature trend of the inverse susceptibility y“^(r). Solid lines: 
a dependence of TciV) oc is expected for P —> 0. For 

strong V, the data are consistent with Tc{V) oc ol{V). Dashed 
line: Tc{V) oc jV'^ represents a good approximation to 
a{V) at 1/ = 8. 


tion strength V as can be inferred from the fact that toa 
is considerably higher than mimp already for 1^ = 3. 

Furthermore, the Curie temperature is seen to decrease 
with increasing V in this regime; Tq drops by about a fac¬ 
tor two when increasing the hybridization strength from 
V = 2 to fo = 3. This can consistently be explained by 
referring to the strong-fo limit where the effective model 
Eq. g applies and where the only energy scale is given 
by the coupling a [see Eq. ([^] which decreases with in¬ 
creasing V. 

This also means that the crossover regime shifts to 
stronger hybridizations strengths with increasing lattice 
dimension or coordination number. This must be kept in 
mind when comparing Tq obtained for different dimen¬ 
sions D at constant V. AtV = 2, the Curie temperature 
does not depend very much on D-. We find Tq = 0.040 
for D = 1, Tc = 0.059 for D = 2 and Tc = 0.064 for 
D = 3. This is easily explained as a balance between 
two counteracting effects, namely an increase of Tc with 
increasing D characteristic for a mean-field theory on the 
one hand and the mentioned shift of the crossover regime 
resulting in a lower Tc on the other hand. 

Obviously, the D = 1 and D = 2 results are not consis¬ 
tent with the Mermin-Wagner theorenP^ which excludes 
spontaneous breaking of the SU(2) spin rotation sym¬ 
metry for H < 2 at finite temperatures. As a matter 
of course, it cannot be satisfied within a static or within 
dynamical mean-field theory since long wave-length mag¬ 
netic excitations do not feed back to the single-particle 
self-energy. We nevertheless expect that the finite Tc 
predicted by DMFT is physically significant even for 


D = 1 (and D = 2) and indicates the onset of ferromag¬ 
netic ordering of the magnetic moments on intermedi¬ 
ate length scales.l^ This corresponds to a thermodynam¬ 
ically stable ferromagnet only if the SU(2) symmetry is 
broken explicitly, e.g., due to the presence of additional 
anisotropic terms in the Hamiltonian. 

For D = 1 and in the low-temperature limit the DMFT 
agrees well with essentially exact data obtained by means 
of the density-matrix renormalization group (DMRG) 
method^ at zero temperature. The extrapolation of the 
DMFT results for the ordered magnetic moments mimp 
and tob to T — 0 perfectly matches with the DMRG 
data, see black squares in Figs. and [T^ As concerns 
the magnetic moment on the A sites, we expect the same 
unconventional T-dependence that has been discussed in 


i.e. 


the context of static mean-field theory in Sec. IID 
an upturn of toa for T —> 0, consistent with the T = 0 
DMRG data, which is induced by the van Hove singu¬ 
larity of the spin-dependent local density of states at the 
Fermi edge. 

For the D = 3 lattice, we have systematically com¬ 
puted Tc as a function of V at fixed U = 8. Results as 
obtained from by linear fits t o th e temperature trend of 
the x“^(T) are shown in Fig. 
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In the weak-coupling limit H —>■ 0, the Curie tem¬ 
perature is expected to be solely determined by the 
strength of the RKKY interaction and thus to scale as 
Tc oc Jrkky oc oc with V. For the strong-coupling 
or HME limit, the only energy scale of the effective low- 
energy theory is given by a and thus Tc (X a (see solid 
lines in Fig. [TI| ). For 1/ = 8, a good approximation is 
Tc oc V~'^, see Eq. ^ and the dashed line in Fig. 
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The Curie temperature is at its maximum Tc,opt ~ 
0.07 for a hybridization strength of about t4pt « 2.5. 
Tc_opt is almost an order of magnitude smaller than the 
maximum Neel temperature of the D = 3 Hubbard model 
at half-filling^ and also an order magnitude smaller than 
typical Curie temperatures of the Hubbard model with 
asymmetric free density of states, as obtained for lower 
fillings by DMFT.I^ The same holds if compared with 
DMFT estimates for the Curie temperature of the stan¬ 
dard periodic Anderson model.l^ 


IV. SUMMARY 

The present study has demonstrated that the 
Anderson-lattice model with a regularly depleted system 
of localized orbitals at every second site supports fer¬ 
romagnetic long-range order which exhibits, depending 
on the hybridization strength U, a high stability against 
thermal fluctuations. The temperature-dependent mag¬ 
netism has been investigated systematically for different 
coupling strengths and electron densities. 

The depleted Anderson lattice model has been consid¬ 
ered beforehand to study f unda mental questions of mag¬ 
netic coupling mechanism^SliO] and to describe artificial 
Hondo systems realized as ultracold atoms trapped in 
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optical lattices.^ It is related to two-dimensional super¬ 
lattices consisting of periodic arrangements of /-electron- 
and non-interacting layer^^ and may be used to describe 
systems of magnetic atoms on non-magnetic metallic sur¬ 
faces where a manipulation of the adatom geometry and 
a precise mapping of magnetic couplings is accessible to 
scanning-tunneling techniques on an atomic scalel^^M^ 

We have employed two different types of mean-field 
approaches: (i) static mean-field theory of the effective 
low-energy model that emerges at strong couplings V 
within fourth-order perturbation theory, and (ii) dynam¬ 
ical mean-field theory of the full model using continuous¬ 
time quantum Monte-Carlo as impurity solver. The 
Curie temperature is obtained by computing the tem¬ 
perature dependence of the magnetic moments as well 
as by the divergence of the homogeneous static magnetic 
susceptibility. The maximal Tc is found at half-filling 
and for intermediate hybridization strengths: 

For weak V, magnetic order is induced by the standard 
effective RKKY interaction between the local magnetic 
moments formed at the correlated impurity sites. For the 
geometry considered and at half-hlling, the RKKY inter¬ 
action is ferromagnetic. The Curie temperature scales 
with in this limit. For strong V, on the other hand, 
the recently proposed inverse indirect magnetic exchange 
also leads to ferromagnetic order. In this limit the im¬ 
purity magnetic moments are Kondo screened and form 
almost local Kondo singlets on a high-energy scale V 
which localize the fraction of conduction electrons not 
taking part in the screening. Those conduction electrons 
develop local magnetic moments which are ferromagnet- 
ically coupled by virtual excitations of the local Kondo 
singlets on an energy scale a [see Eq. (|^]. Therefore, Tq 
scales with a ^ V~^ for fixed U in this limit. 

While the numerical data obtained for different V ap¬ 
pear to be consistent with the expected trends, it turned 
out to be very difficult to reach the extreme limits K —>■ 0 
and K —?► oo characterized by pure RKKY or IIME cou¬ 
pling, respectively, as the energy scale given by Tc be¬ 
comes too small. As concerns the strong-coupling limit, 
we conclude that a perfect linear scaling of Tq with a can 
only be expected for still stronger hybridization strengths 
V that are not accessible to DMFT with the presently 
used impurity solver. This also implies that a direct com¬ 
parison of the DMFT results for Tq with those obtained 
by static mean-field theory applied to the effective low- 
energy model is not meaningful. From Fig.j^we can infer 
that the latter would predict a Curie temperature which 
is by two orders of magnitude higher than the DMFT 
result for K = 7 in Fig. El This could indicate that the 
strong-coupling limit is not yet reached but could also be 
ascribed to strong local fluctuations reducing Tq which 
are accounted for within the dynamical but not in the 
static mean-field theory.l^ 

On the other hand, the maximum Tq opt ~ 0.07 found 
for intermediate V is well accessible to DMFT and sur¬ 
prisingly high, in view of the fact that the magnetic cou¬ 
pling is mediated indirectly only. Compared to DMFT 


estimated^^ll^ of critical (Neel or Curie) temperatures 
in the Hubbard or periodic Anderson model with a 
dense system of correlated impurities, it is about an or¬ 
der of magnitude lower. The optimal intermediate hy¬ 
bridization strength where Tq is at its maximum is given 
by I4pt ~ 2.5, i.e., clearly stronger than the nearest- 
neighbor hopping t = 1. 

One should note that DMFT applied the depleted An¬ 
derson lattice can be expected to be much more reliable 
than for the dense case. In fact, perfect agreement with 
numerically exact DMRG data is observed in the low- 
temperature limit. We are therefore convinced that this 
study provides quantitative insight into the physics and 
contributes to the fundamental understanding of mag¬ 
netic order of correlated orbitals coupled indirectly by 
conduction electrons. 
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Appendix A: Low-temperature behavior of the 
magnetization at half-filing 


Here we derive the low-temperature behavior of 
dmx/dT at half-filling, i.e., Eqs. (16) and ( [T7| . Due to 
particle-hole symmetry at half-hlling (^ = 0), we have 




(Al) 


and 


a 


r?t(k) = -r?^(k) = --(D-Q;)7^(k) (A2) 


Using Eq. (15), we immediately have 

a(D-Q;(T)) 72 (k) 


1 


ruA 


= — > tanh • 

La ^ 


fceBZ^ 


4T 


(A3) 


In the thermodynamical limit, T i—>■ oo, the k-sum can 
be replaced by one-dimensional integration. 


ruA 



dujp^^\uj) tanh 


ajP - Q^{T))uj^ 

4T 


(A4) 


with the weight function 

p(D) = E ^ S{uj - 7(k)). (A5) 


We note that and p^^^(a;) are Hnite at the Eermi 

edge (at co = 0), while ^ ln|a;| diverges. Since 
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= p^^\uj) and since Qi{T) ~ 0 for low tem¬ 
peratures, Eq. (|A4|) implies 


dniA 




UJ 


cosh^ ■ 


After changing the integration variable, we obtain 
dmA _ 1 f°° _(D) ( r, I Tx \ ^/x 


dT 


.if-'-fA) 


cosh"^ 


(A6) 


(A7) 


for low T. Since p^^\uj) is regular at w = 0 for _D = 1 
and D = 3, we have 


diTiA 


T-J-O 1 

r=o ^/aDT 


(£) = !,£) = 3) (A8) 


and thus Eq. (16). Eurthermore, with p^'^\yj2T/a) 
ln(2T/a) we get 


dm A T-s-o ln(r/a) 


dT 


y/aDT 


{D = 2) 


(A9) 


i.e., Eq. ( [T7| . The low-temperature behavior of mA is 
thus governed by the weight function p^^'>(uj) at low w. 
Analogously, one may also relate the low-temperature be¬ 
havior of toa to the van Hove singularity of the spin- 
dependent tight-binding density of states corresponding 
to the mean-field dispersion Eq. (13). 
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